Influence of entrance channels on formation of 
superheavy nuclei in massive fusion reactions 
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Abstract 

-t-p* , Within the framework of the dinuclear system (DNS) model, the production cross sections of 

O I superheavy nuclei Hs (Z=108) and Z=112 combined with different reaction systems are analyzed 

, systematically. It is found that the mass asymmetries and the reaction Q values of the combina- 

I ■ tions play a very important role on the formation cross sections of the evaporation residues. Both 

> ; 

■rj" I methods by solving the master equations along the mass asymmetry degree of freedom (ID) and 
G^ . 

0^ ■ along the proton and the neutron degrees of freedom (2D) are compared each other and with the 
Tj- ■ available experimental results. 
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^ ' The synthesis of heavy or superheavy nuclei (SHN) is a very important subject in nuclear 
physics motivated with respect to the island of stability which is predicted theoretically, and has 
obtained much experimental research with fusion-evaporation reactions [H Combinations with 
a doubly magic nucleus or nearly magic nucleus are usually chosen owing to the larger reaction 
Q values. Six new elements with Z=107-112 were synthesized in cold fusion reactions for the 
first time and investigated at GSI (Darmstadt, Germany) with the heavy-ion accelerator UNILAC 
and the SHIP separator [H [3]. Recently, experiments on the synthesis of element 113 in the 
'^''Zn+^'^^Bi reaction have been performed successfully at RIKEN (Tokyo, Japan) [Ij. However, it 
is difficulty to produce heavier SHN in the cold fusion reactions because of the smaller production 
cross sections that are lower than 1 pb for Z > 113. The superheavy elements Z=113-116, 118 
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were synthesized at FLNR in Dubna (Russia) with the double magic nucleus "^^Ca bombarding 
actinide nuclei O El [7]. New heavy isotopes ^^^Db and ^^^Bh have also been synthesized at 
HIRFL in Lanzhou (China) [8]. Further experimental works are necessary in order to testify the 
new synthesized SHN. A better understanding of the formation of SHN in the massive fusion 
reactions is still a challenge for theory. 

In this letter, we focus on the influence of the entrance mass asymmetry and the reaction Q 
value of projectile-target combinations on the production cross sections of superheavy residues. 
In the DNS model, the evaporation residue cross section is expressed as a sum over partial waves 
with angular momentum J at the centre-of-mass energy Ec.m. O [IHl [TT] . 

Jmax 

(TEBiE^.m.) = ^ E (2-/ + l)2^(^c.n^., J)Pcn{E,,^,, J)Wsur{E,,^,, J). (1) 

Here, T{Ec,m., J) is the transmission probability of the two colliding nuclei overcoming the Coulomb 
potential barrier in the entrance channel to form the DNS. The Pcn is the probability that the 
system will evolve from a touching configuration into the compound nucleus in competition with 
quasi-fission of the DNS and fission of the heavy fragment. The last term is the survival proba- 
bility of the formed compound nucleus, which can be estimated with the statistical evaporation 
model by considering the competition between neutron evaporation and fission We take the 
maximal angular momentum as Jmax = 30 since the fission barrier of the heavy nucleus disappears 
at high spin |12j . 

In order to describe the fusion dynamics as a diffusion process along proton and neutron degrees 
of freedom, the fusion probability is obtained by solving a set of master equations numerically 
in the potential energy surface of the DNS. The time evolution of the distribution probability 
function P{Zi, Ni, Ei,t) for fragment 1 with proton number Zi and neutron number A^i with 
excitation energy Ei is described by the following master equations |13j . 

dPiZ,,N,,E,,t) 



Y.Wz„N,;Z[,NAt) [dz„N,PiZ[,N,,E[,t) - dz'^,N,PiZi,NuEut) 
zi 



+ 



T.Wz^,Nr,ZuNi{t) [dz„N,P{ZuKE[,t)-dz„N[PiZi,Ni,Ei,t) - 

K 

A«^(e(t)) + A^*^(e(t))] P(Zi, N,,E^,t). (2) 

Here Wzi,Ni]Z[,Ni iWzi,Ni;Zi,N[) is the mean transition probability from the channel {Zi, Ni, Ei) 
to {Z[, Ni, E[) (or {Zi, Ni, El) to {Zi, N[, E[)) , and dz^,Ni denotes the microscopic dimension 
corresponding to the macroscopic state {Zi,Ni,Ei). The sum is taken over all possible proton 
and neutron numbers that fragment Z[, N[ may take, but only one nucleon transfer is considered 



in the model with Z[ = Zi ± 1 and = A^^i ± 1. The excitation energy Ei is determined by 
the dissipation energy from the relative motion and the potential energy surface of the DNS. 
The motion of nucleons in the interacting potential is governed by the single-particle Hamiltonian 
[H [ID]- The evolution of the DNS along the variable R leads to the quasi-fission of the DNS. The 
quasi-fission rate A''-^ and the fission rate A-^** can be estimated with the one-dimensional Kramers 
formula [H [11]. 

In the relaxation process of the relative motion, the DNS will be excited by the dissipation of 
the relative kinetic energy. The local excitation energy is determined by the excitation energy of 
the composite system and the potential energy surface of the DNS. The potential energy surface 
(PES) of the DNS is given by 



UiZi,Ni,Z2,N2;J,R;Pi,P2,0u92) = B{Zi,Ni) + BiZ2,N2) - B{Z,N) + vZt {J) 



+ 



A^i, Z2, A^2; J, R; A, /?2, ^1, 92) (3) 

with Zi + Z2 = Z and A^^ + A^2 = A^- Here B{Zi,Ni){i = 1,2) and B{Z,N) are the negative 
binding energies of the fragment (Zj, Aj) and the compound nucleus {Z, N), respectively, in which 
the shell and the pairing corrections are included reasonably. The V^^^ is the rotation energy 
of the compound nucleus. The /3j represent the quadrupole deformations of the two fragments. 
The 6i denote the angles between the collision orientations and the symmetry axes of deformed 
nuclei. The interaction potential between fragment {Zi,Ni) and {Z21N2) includes the nuclear. 
Coulomb and centrifugal parts, the details are given in Ref. [TU]. In the calculation, the distance 
R between the centers of the two fragments is chosen to be the value which gives the minimum of 
the interaction potential, in which the DNS is considered to be formed. So the PES depends on 
the proton and neutron numbers of the fragment. In Fig.l we give the potential energy surface 
in the reaction ^'^Si+^^^Cf as functions of the protons and neutrons of the fragments in the left 
panel. The incident point is shown by the solid circle and the minimum way in the PES is added 
by the thick line. The driving potential as a function of the mass asymmetry that was calculated 
in Ref. [HI HO] is also given in the right panel and compared with the minimum way in the left 
panel. The driving potential at the incident point in ID PES is located at the maximum value, so 
there is no the inner fusion barrier for the system, which results in a too large fusion probability. 
Therefore, we solve the master equations within the 2D PES to get the fusion probability for the 
systems with larger mass asymmetries. 

The formation probability of the compound nucleus at the Coulomb barrier B (here a barrier 



9. 



distribution f{B) is considered) and for angular momentum J is given byjOl ITU] 

Zbg ^bg 

Pcn{E,,^.J.B)= J2 E PiZ,,N,,E,,T,^t{E,,^,,J,B)). (4) 

Zi=l Ni=l 

We obtain the fusion probability as 

PcNiE,.m., J) = J f{B)PcN{E,,^. J, B)dB, (5) 

where the barrier distribution function is taken in asymmetric Gaussian form. 

The survival probability of the excited compound nucleus cooled by the neutron evaporation 
in competition with fission is expressed as follows: 

"^•-(^-- ^) = '^^'c' ^) n ( r„(E.-'j5+'r'U.j) ). ' 

where the -E'^Af' excitation energy and the spin of the compound nucleus, respectively. 

The E* is the excitation energy before evaporating the ith neutron, which has the relation 

E*^, = E* - B: - 2T,, (7) 

with the initial condition E^ = E^j^. The energy B^ is the separation energy of the ith neu- 
tron. The nuclear temperature Tj is given as E* = aTf — Ti with the level density parameter a. 
P{Eq]^,x, J) is the realization probability of emitting x neutrons. The widths of neutron evap- 
oration and fission are calculated using the statistical model. The details can be found in Refs. 

mm 

With this procedure introduced above, we calculated the evaporation residue excitation func- 
tions using the ID and 2D master equations in the reaction ^^Ca+^^^U as shown in Fig.2 rep- 
resented by dashed and solid lines, respectively, and compared them with the experimental data 
performed in Dubna [H] and at GSI [15]. The GSI results show that the formation cross sections 
in the 3n channel at the same excitation energy with 35 MeV have a slight decrease, which are in 
a good agreement with our ID calculations. In the whole range, the 2D calculations give smaller 
cross sections than ID master equations owing to the decrease of the fusion probability. For the 
considered system, the value of the PES at the incident point is located at the line of the min- 
imum way. So the ID master equations can give reasonable results. However, for the systems 
with larger mass asymmetries and larger quadrupole deformation parameters, e.g. 
^^Ne-|-^^^Pu, etc, the ID master equations give too large fusion probabilities. 

The synthesis of heavy or superheavy nuclei through fusing two stable nuclei is inhibited by 
the so-called quasi-fission process. The entrance channel combinations of projectile and target 
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will influence the fusion dynamics. The suppression of the evaporation residue cross sections for 
less flssile compound systems such as ^^^Ra and ^so^jj^ when reactions are involved in projectiles 
heavier than ^^C and ^^O was observed in Refs. [16]. The wider width of the mass distributions for 
the flssion-like fragments was also reported in Ref. [T7|. In Fig.3 we calculated the transmission 
and fusion probabilities using the 2D master equations for the reactions ^^S+^'^^U, ^^Fe+^''^Pb and 
^^^Xe+^^^Xe which lead to the same compound nucleus ^^^Hs formation. The larger transmission 
probabilities were found in the reactions ^^Fe+^''^Pb and ^^^Xe+^^^Xe owing to the larger Q values 
(absolute values). Smaller mass asymmetries of the two systems result in a decrease of the fusion 
probabilities. The evaporation residue excitation functions in ln-5n channels are shown in Fig. 4. 
The competition of the capture and the fusion process of the three systems leads to different trends 
of the evaporation channels. The 3n and 4n channels in the reaction ^^S+^^^U, In and 2n channels 
in the reaction ^^Fe+^°'^Pb are favorable to produce the isotopes 269,268jjg 27i,270jjg_ Although 
the system ^^^Xe+^'^^Xe consists of two magic nuclei, the higher inner fusion barrier decreases the 
fusion probabilities and enhances the quasi-flssion rate of the DNS, hence leads to the smaller cross 
sections of the Hs isotopes. The upper limit cross sections for evaporation residues cr(i_3)„ <4 pb 
were observed in a recent experiment [18], which are much lower than the ones predicted by the 
fusion by diffusion model [19]. In the DNS model, the larger mass asymmetry favors the nucleon 
transfer from the light projectile to heavy target, and therefore enhances the fusion probability of 
two colliding nuclei. 

The superheavy element Z=112 was synthesized at GSI with the new isotope '^'^^112 in cold fu- 
sion reaction ™Zn+^°^Pb [20] and also fabricated with more neutron-rich isotopes 282,283 ]^]^2 in ^^Ca 
induced reaction ^^Ca+^^^U. We analyzed the combinations ^^Si+^^^Cf, ^^S+^^^Cm, '^OAr+^^^'Pu 
and ^^Ca+^^^U which lead to the production of new isotopes of the element Z=112 between the 
cold fusion reactions and the ^^Ca induced reactions as shown in Fig.5. The 2n, 3n and 4n channels 
in the reaction ^°Si+^^^Cf, and the 4n channel in the reaction ^^S-|-^^°Cm have larger cross section 
to produce new isotopes due to the larger fusion probabilities of the two colliding nuclei. 

In summary, we systematically analyzed the entrance channel effects of synthesizing SHN using 
the DNS model. The systems with larger entrance mass asymmetry and larger reaction Q value 
can enhance the capture and fusion probabilities of two colliding nuclei. Calculations were carried 
out for the reactions ^^S+^'^^U, ^^Fe-|-^°^Pb and ^^^Xe+^^^Xe which lead to the same compound 
nucleus formation. The 2n, 3n and 4n channels in the reaction ^''Si+^^^Cf, and the 4n channel 
in the reaction ^^S+^^°Cm are favorable to synthesize new isotopes of the element Z=112 at the 
stated excitation energies. 
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Figure 1: The potential energy surface of the DNS in the reaction ^°Si+^^^Cf as functions of the 
protons and neutrons of the fragments (left panel) and the mass asymmetry coordinate (right 
panel). 




Figure 2: Comparison of the calculated evaporation residue excitation functions using the ID and 
2D master equations with the available experimental data in the reaction ^^Ca+^'^^U. 
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Figure 3: Calculated transmission and fusion probabilities as functions of the excitation energies 
of the compound nucleus for the reactions ^^S+^^^U, ®^Fe+^°^Pb and ^^^Xe+^^^Xe. 




Figure 4: Comparison of the calculated evaporation residue cross sections in ln-5n channels using 
the 2D master equations for the reactions ^^S+^^^U, ^^Fe+^^^Pb and ^^^Xe+^^^Xe. 




Figure 5: The same as in Fig.4, but for the reactions ^"Si+^^^Cf, ^^S+^^^Cm, ^OAr+^^'^Pu and 
^^Ca+^^^U leading to the formation of the element Z=112. 
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